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We outline a scheme to accomplish measurements of a solid state double well system (DWS) with 
both one and two electrons in non- localised bases. We show that, for a single particle, measuring the 
local charge distribution at the midpoint of a DWS using an SET as a sensitive electrometer amounts 
to performing a projective measurement in the parity (symmetric/antisymmetric) eigenbasis. For 
two-electrons in a DWS, a similar configuration of SET results in close-to-projective measurement 
in the singlet /triplet basis. We analyse the sensitivity of the scheme to asymmetry in the SET 
position for some experimentally relevant parameter, and show that it is realisable in experiment. 
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I. INTRODUCTION 

In this paper, we present a scheme for performing measurements of one- and two-electron double well systems 
(DWSs) in non-local bases. The principle idea behind both schemes is to detect the small charge differences at the 
midpoint of the DWS for different electronic states using a single electron transistor (SET) as a sensitive electrometer. 
set's have been demonstrated to be highly sensitive electrometers, with sensitivities of a few iie/VHz We 
therefore expect that they could detect fluctuations of around 1% of an electron charge in around 10-100 ns. 

For the one electron case, the measurement is in the parity eigenbasis. This is interesting since it results in projection 
onto non-local states. It has been shown theoretically and experimentally that decoherence is slower for an 
evenly biased DWS, so this kind of measurement may be more robust against decoherence. 

For a two electron system we show that this measurement approximately projects the DWS onto the singlet (even) 
and triplet (odd) subspaces. It is therefore a method for performing spin sensitive detection using electrometers, 
which is important for readout of certain quantum information processing schemes 0. 

The paper begins in section^ with a short, generic discussion of the microscopic model of a few electron system 
interacting with an idealised SET, itself in contact with a continuum of lead modes. Following this, section ITTll deals 
with the single electron case, and section IIVI deals with the two electron case. Within each of these two sections 
we develop the measurement Hamiltonian from microscopic considerations, from which we derive measurement and 
mixing times. After showing in each case that the measurements work in principle, we estimate the effects of a 
significant problem in the fabrication of this device, namely the precision with which the SET island must be placed 
in the midplane of the DWS. Section El concludes the paper with a discussion about experimental implementation of 
the scheme. 

►SET island 
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FIG. 1: Schematic of physical system under consideration 
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II. SYSTEM 



We consider the device pictured in Fig. ^ consisting of a SET placed in the mid-plane of the DWS, in order 
to be sensitive to the charge at the midpoint. This charge differs between symmetric and antisymmetric spatial 
wavefunctions, and we analyse a scheme to measure this difference in order to effect projective measurements onto 
the parity eigenspaces. 

For the purposes of this paper, we assume the SET island has only a single accessible energy level, which is reasonable 
if the island is small and the difference between Fermi energies in the leads is less than the charging energy of the 
island. We model the DWS with a Hubbard Hamiltonian, with only a single spatial wavefunction per well, \L) and 
\B) for the left and right wells. We assume that the Hilbert space for the system is therefore two-dimensional, which 
is reasonable if higher excited state are inaccessible due to the strong confinement of the quantum dot potentials. 
Therefore, each well may be populated by at most two electrons, in different spin configurations. 

The Hamiltonian for a system of interacting electrons is given by 



ij.17 ijlm^aa' 

where Ci^a is a (fermionic) annihilation operator for an electron in spatial mode i and spin a G {t, i} Q| and 

H,j = j d?v<t>:{r){-feV^/2m + U{r))<t>j{r) (2) 
dV |d3r2</)*(ri)0,(ri)F(|ri-r2|)C(r2)0™(r2), (3) 



where (f>i{v) = {y\c\^\) is the spatial wavefunction for mode i (assuming both spin states have the same spatial wave- 
function), I) is the 'vacuum' state, with no quasi-particlc excitations, m is the effective mass, U{r) is the electrostatic 
confining potential and V{r) = ~ f ^^'-^ Coulomb potential. 

We choose a basis set such that Hij is diagonal, which we truncate to the lowest two eigenstates for the DWS, 
and a single state, s, on the SET. In the absence of an external bias between the wells, this corresponds to taking 
I, j S {-|-, — , s}, where ± are the symmetric and antisymmetric superpositions of the localised, single-particle states, 
referred to here as the parity eigenbasis. That is |±) = cl|) ~ \'^)MR) ^\ic first term of Eq. dTJ becomes 

^ I / ±1/ ^2(1±{/?|L>) ^ *-* 

-^(n+o- — n„(j) + oJofis, where hia = c\^Cia and A = — H is the tunnelling rate of a localised electron, 



which we estimate in section Ivl 
Expanding Eq. gives 



-ffTot = iJisland + ^^Icads + -H^tun + i^Hub + V , (4) 



where i/nub is the Hubbard Hamiltonian for the DWS, V is the Hamiltonian for the interaction between the DWS 
and the SET and 

-ff island = ^O^s, (5) 

-ffleads = ^ + firfc), (6) 

k 

Htnn = Y^Tikc\^.Cs +Trkclf.Cs + W.C., (7) 
k 

Here ljq is the island energy level, in the absence of interactions with the double well potential, lJ]^ are the energies of 
densely spaced lead modes, / and r denote the left and right leads respectively, Tn^r)k are the corresponding tunnelling 
rates between mode k in lead l{r) and the SET island. 

To compute V we assume that the wavefunction for electrons on the SET island vanishes in the region where the 
wavefunction for the electron on the DWS has support, and vice versa, so that for instance (^!^(r)(/)s(r) = 0. This 
assumption is a good one for systems such as the Kane proposal Q, or Na+ in Si [3, where the tunnelling rates 
between SET island and the DWS are negligible. 

The consequence of this assumption is that if any index in Vijmn is s, then Vijmn is zero unless i = j = s or 
m = n = s, where s labels an electron on the SET island. Therefore the only Coulomb terms that contribute to the 
interaction between the SET island and the DWS are given by 

V = ^ Vssijflscl^Cja =hs<E) -ffmoas, (8) 
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where we have ignored the spin degree of freedom on the SET island. There arc four distinct terms of this form 
that need to be included, for i,j taking the four possible assignments of + or — . If the physical arrangement 
of double well and island as shown in Fig. ^ is symmetric about a line bisecting the double well potential, then 

e = Vss-i = Vss h = 0, as discussed in appendix^ Asymmetry results in non-zero e, and for the development of 

this section, we assume that it is symmetric, so that e = 0. Thus, to describe a DWS interacting with an SET we 
need to specify i?Hub and Hmcas- 

III. PARITY MEASUREMENT FOR SINGLY OCCUPIED DWS 

The system we consider in this section consists of a single electron shared between two wells, so we ignore spin 
indices. We now establish the feasibility of performing a measurement in the parity eigenbasis of a single electron 
shared between the two wells. 



A. Derivation of measurement Hamiltonian 

For a single electron on an unbiased double well, the interaction terms of Eq. vanishes, so the Hubbard Hamil- 
tonian for the single particle system is given by 

Hnnh = An+. (9) 

We have used the fact that, for a single-particle, two-level system fi^ + fi- = I is the identity, and have also discarded 
terms proportional to /. In an alternate notation, we identify n+ with the ax = + \R){L\ operator. 

For the symmetric case, 

V ^hs® {Vss++n+ + Vss — = Shsh+, (10) 

where S = Vss++ — Vss and we have again used the single particle, two-level system identity 77.+ + fi- = /. It 

is evident from Eqs. ^ and (|10|l that the measurement Hamiltonian and the system Hamiltonian have the same 
energy eigenstates. Therefore, the measurement process will non-destructively project onto energy eigenstates of 
the DWS, which are the delocalised symmetric and antisymmetric wavefunctions. That is, it is a QND (quantum 
non-demolition) measurement, which simplifies the analysis greatly. 

B. Master equation for symmetric system 

As discussed in appendix^ the master equation for the DWS and SET island is given by 

R{t) = - t[HHuh + ^^island, R{t)] + {jl + JrMh^cl]R{t) + j',V[h+Cs]R{t) + 7,'P[r7+ct ] i?(t) , (11) 

where R is the density matrix for the DWS and SET island a.ndV[A]B = J[A]B-A[A]B = ABA^-^iA'^AB+BA^A). 
We assume that the reduced state of the SET island is diagonal in the number representation, so 

i^^Po® |0).(0|+pi® |1),(1|, (12) 

The reduced density matrix of the DWS is given by p{t) = TTs{R{t)} — po{t) + pi{t). In appendix IXlwe solve Eq. Hll() 
and the steady-state reduced density matrix of the DWS is 

Tr4i?(oo)} = p{^) = (p++(0) + pt+mh+ + {p^-{0) + p^-mn-, (13) 

where p'"^ = {p\p\q). This shows that the diagonal elements of the system density matrix are unchanged, whilst the 
off-diagonal elements have completely decayed, which corresponds to QND measurement. 

C. Measurement time 

Since the measurement for the single electron case is a QND measurement, we can straightforwardly calculate 
the expected currents and measurement times for the device. If we choose fi{i^o) = fii^^Jo + S) = I ~ fr{^a) and 
fr{(jJo + S) = 0, as shown in Fig. 12 then the current through the SET is sensitive to the state of the DWS. In 
this configuration, the SET island is conducting when the DWS is in a symmetric state, so a current 7+ flows, and 
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FIG. 2: SET island energies relative to lead Fermi level, depending on state of DWS. 



FIG. 3: Geometry of asymmetric system. 

non-conducting when in an antisymmetric state, so i_ = 0. Thus the measurement amounts to distinguishing the 
currents i+ and i_ = through the SET. 

For the symmetric configuration, the rates at which electrons hop on and off the SET island, are given by 7; 
and 7r. Thus the rate of transport of electrons through the SET island is 1/(7^"^ + 7^^). The current is therefore 

i+ = e/(7r^ +77^)- 

The measurement time is then the time required to distinguish two currents, i+ from i_ = in the presence of 
shot noise. Since the transport of electrons through the SET is Poissonian, the probability of detecting zero electrons 
tunnelling in a time r through the SET island, given that a mean current i+ is flowing is given by Po(T) = I ^'^i ^ . 
We therefore conclude that the measurement time is tmeas ~ 7;~^ + 7^^! since the probability of not detecting a 
tunnelling event in (a few multiples of) this time is small. This agrees with the measurement determined from the 
decay rate of off diagonal elements of the system density matrix, found in appendix fXl 



D. Effects of asymmetry 

We now turn to the important practical issue of how precisely the SET island needs to be placed with respect to 
the centre of the double well potential. Here we assume there is some asymmetry, which may, for example, arise from 

fabrication, so that e — I^sh 7^ 0. We will estimate the magnitude of this quantity later, but first we will determine 

the effect of the extra terms in the Hamiltonian that arise. Including this term in V gives 

V" ® + e(c^c_ + cLc+)^ . (14) 

This shows that the asymmetry rotates the measurement basis by an angle = tan~^(e/(5) away from the parity 
eigenbasis, that is, the preferred basis for the measurement is {cos((/))|-|-) + sin((/))|— ), — sin((?!))|+) + cos((/i)|— )}. We 
therefore require that e ^ 5 in order that the asymmetry have negligible effect. 
In Appendix IbI we estimate e and 8 to be 

' |rL - r^-^P ^ |rL - r7"T ^""^ ^ ^ '^^^'^^ " Vl - rl^'"| ) ' ^^^^ 

where tl is the 'centre-of-mass' of the left well, = {xs, ?/s, Zg) is the position vector of the SET island and we choose 
the origin to be at the midpoint of the DWS. Note that = for a symmetrically placed island, and is assumed to 
be small. 

The condition that e ^ is therefore satisfied if | w {l\R) |ri,-rj^"'p ^ This is a tight constraint, since it 
requires that the asymmetry, quantified by ^^_^5ym|2 , be much less than the overlap of the localised wavefunctions 
(L|i?). With the help of a J-gate (as referred to in Kane Q) (-^1^) may be made as high as 0.03 0, and assuming 
a typical scale of device of jr^^ — r^j 30 nm, the elements of the SET and DWS would likely need to be made with 
a precision of 1 nm or less, which seems difficult with current technology. 
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This issue may not be so significant for electrostatically defined dots, since the position of the SET island and DWS 
may be changed by variation of surface gate potentials. It is a serious problem for donor systems with SET's grown 
by metallic deposition, since the location of the donor atoms and SET island are fixed during fabrication. 

E. Mixing time 

Asymmetry in the placement of the SET island induces mixing in the state of the DWS, so there is a mixing time 
associated with asymmetry. The calculation of the mixing time is somewhat lengthy, but not difhcult. We derive an 
unconditional master equation for the density matrix of the DWS and SET, R. The solution to the master equation 
has exponentially decaying terms, with different time constants. For the sake of simplicity, here we present the results 
of the calculation, and leave the details to appendix^ which follows from the results of appendix|0 Taking 7,; ~ 7, 
the most rapidly decaying term gives the measurement time, tmeas = 2/7, which is unchanged from the symmetric 

case (to within 0(e^)). The slowest decaying term gives the measurement induced mixing time, ~ ^'^^2^'' j^i so 
the condition for a good measurement is that t^^i^ 3> tmcasj which occurs if (5 + A ^ e, in agreement with above. 
When e = 0, the measurement is QND, imix = 00. 

When e 7^ 0, the dynamics arc divided into two regimes, t <^ imix and t ^ tmix- In the energy eigenbasis, for short 
times, t ^ <mix, the diagonal elements of the density matrix are almost unaffected whilst the off-diagonal elements 
decay at a rate 7/2. This corresponds to the process of projecting the DWS onto its energy eigenstates. Over longer 
times, t ~ imix, the diagonal matrix elements decay to their steady-state values, given by Eq. ljD6|l . corresponding to 
almost complete relaxation of the DWS to its ground state |— ). 

F. Summary 

Performing a 'non-local' measurement by projecting onto the parity eigenbasis of a DWS is in principle possible. 
An SET placed in the midplane of the DWS is sensitive to the local charge at the midpoint of the system, which 
depends on the parity of the state. 

Estimates of the sensitivity of such a device to asymmetry in the placement of the SET island suggest that it would 
need to be placed with a precision better than {L\R) sa 3% of the typical device dimensions. For a donor based 
system, this is around 1 nm, which has been demonstrated recently for P+ donors in Si 

The scheme may still be usefully applied in the case where e ^ where the preferred measurement basis is close 
to the parity eigenbasis. Applying a bias across the DWS, so that one well is at a higher potential than the other, 
would allow one to rotate the double well energy eigenbasis onto the measurement eigenbasis. 

The precise placement of the SET island in an electrostatically defined system may not present serious problems, 
since it may be moved about after fabrication by varying surface gate voltages. In section we give estimates of 
expirementally accessible parameter values, and show that this proposal is experimentally viable. 

IV. SINGLET/TRIPLET MEASUREMENT FOR DOUBLY OCCUPIED DWS 

We consider here a DWS populated with two electrons, shown schematically in Fig. ^ As in the single electron 
parity measurement scheme, the SET island is placed in the midplane of the DWS, so that it is sensitive to the charge 
distribution on the DWS. 

The physical principle that we exploit is the fact that a pair of electrons in a triplet state are Pauli blocked from 
being simultaneously at the origin of the DWS, so the probability amplitude to find two electrons in a triplet at the 
origin is zero. In contrast, this restriction does not apply to a pair in a singlet state, so there is a non-zero probability 
amplitude to find two electrons in a singlet near the origin. Thus there is a small variation in the local charge density 
at the origin between singlet and triplet states, which can in principle be measured to distinguish these subspaces. 
This has some similarities to another singlet-triplet measurement scheme |l2j . 

Such systems may be used to implement quantum information processing tasks. In certain instances it is important 
to distinguish whether the two electrons are in a singlet state or a triplet state, thereby providing information about 
their spin state, e.g. distinguishing a state from the triplet states is necessary in the three-in-one encoding scheme 
developed by DiVincenzo et al. 0. 



6 



A. Derivation of measurement Hamiltonian 

In this system, there is one singly occupied singlet state and three singly occupied triplet states, which are given by 

(16) 

(17) 
(18) 

(19) 

The spatial wavcfunction of the singlet state is clearly symmetric, whilst the spin wavefunction is antisymmetric. The 
converse is true for the triplet states. Written in the parity eigenbasis, it is clear that there are different charge densities 
between the wells depending on the state: the singlet state is a superposition of terms with non-zero amplitude to 
find either zero electrons (both in the |— ) state, with zero charge density at the midpoint) or two electrons (both in 
the 1+) state, with non-zero charge density at the midpoint) to exist between the wells, whilst the triplet states have 
an amplitude to find only a single electron (only one electron in the |+)) state to be located at the midpoint. There 
are also two doubly occupied states given by 



\S) 






























,)l>. 



\D^ 



T"-+i 



(20) 
(21) 



As discussed in sectional we need to specify the DWS dynamics, given by i?Hub, as well as the interaction between 
the DWS and the SET, given by iJmcas- In the parity eigenbasis, iJnub is given by 



Huuh 



A/2(n+| + h+i - h-1 
+ U/2 (^(n+i + n_|)(n+^ + n_^) + (ct_^c_| + cl^c+^){cl^c^i + c^^c+j^)) , 



(22) 



where A = — H as defined earlier and U = Vllll is the double occupation Coulomb energy. The triplet 

states are eigenstates of the two-site Hubbard Hamiltonian, so decouple from the other states. With respect to the 
ordered sub-basis {15), l^*-)} the matrix for the Hubbard Hamiltonian is 



Huuh — 



-U A 
A 




(23) 



We now turn to the interaction between the two electron system and the nearby SET island populated with at 
most one electron. There are four distinct terms in Eq. ||HJ) that need to be computed, for i,j taking the four possible 

assignments of -I- or — . The triplet states are once again eigenstates of i?meas, with eigenvalue Vss++ + V^s Again 

with respect to the ordered sub-basis {\S), |-D+), |-D_)}, the matrix representation of -ffmcas is 



H„ 



a 6 
S a 2e 
2e a 



(24) 



where a = Vss++ + Vss , S = Kss++ — Kss-- and e = T4s-i as defined earlier. As discussed in Section llll Dl if the 

physical arrangement of DWS and SET island as shown in Fig. ^ is symmetric about a line bisecting the double well 
potential, then e = 0. Asymmetry results in e 7^ 0. 

Finally, we include the internal dynamics of the leads and of the island, which is assumed to have at most a single 
electron, as well as tunnelling between the leads and the island. The Hamiltonians for these parts of the complete 
system are given in Eqs. (0) to (0). We now have the ingredients for the model Hamiltonian of the DWS, SET island 
and leads 



J^Tot = Huuh + V + His 



land 



(25) 
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B. Measurement of symmetric configuration 

111 this section, we will assume that the SET is placed symmetrically with respect to the DWS, so e = 0. The 
triplet states and \D_) are energy eigenstates of both -ffnub and -ffmcas- Similarly, since 5 and A are small, IS") is 
approximately an eigenstates of both i^Hub and -ffmeas- The induced charge on the SET island is different for the 
triplet states compared with the singlet, due to their slightly different charge configurations, so different currents flow 
through the SET depending on the subspace the electron pair is in. Thus, distinguishing distinct currents through the 
SET yields (approximately) QND projective measurements onto the singlet or triplet subspaces. This is analogous to 
the situation described for the single particle in section HTll for doing QND measurements in the parity eigenbasis. 

The induced shift in the SET island energy depends on the state of the DWS. We can calculate the SET island 
energy shift by imagining the DWS in a given state, then adiabatically turning on the SET-DWS interaction, V . 
Physically, this corresponds to slowly bringing the occupied SET island close to the DWS, and observing the change 
in energy of the total system during this process. Comparing the adiabatic energy shift for a DWS in a singlet with 
the energy shift for a DWS in a triplet gives the differential shift of the SET island between the singlet and a triplet 
states. 

A triplet state is an energy eigenstates of i/nub + aiJmeas, where < a < 1 is the adiabatic parameter controlling 
the coupling strength. Therefore the adiabatic variation of the coupling does not change the eigenstate, just the 
eigenenergy, 

AEt,, = {T^MHnnb + Hmcas)|Tn) - (Ttt K^^Hub) [Ttt) = a. (26) 

For a singlet state, which is almost, but not quite, an eigenstate of -ffuub + Qfi/mcas, we can estimate the induced 
shift by calculating the same adiabatic energy shift for the ground state \Sa) = l-S*) — ^^^|Z3+) which is very close 
to the singlet. We find 

AEs ~ (5i|(HHub +i?mcas)|^l) - (^0 | i^Hub | ^o) = a - S{6 + 2A) /U. (27) 

Thus, the difference between the SET island energy for the triplet and singlet state is AE = AEt — AEs = 
S{S + 2A)/U, which corresponds to a differential induced SET island energy depending on the state of the DWS. It 
should therefore be possible to arrange the lead energies so that the SET current also depends on the state of the 
DWS. By tuning the lead chemical potentials so that /i/ > Et > Hr > Es, current flows through the SET if the 



DWS is in the triplet subspace, but does not for the subspace {\S), jZ?-!-)}. This configuration is shown in Fig. 5(a) 



where Et = ^^^6.7,8 lies in between the lead Fermi energies, and Es = uj2 lies below the Fermi energies. The other 



levels shown in Fig. 5(a) represent possible inelastic transitions as lead electrons tunnel onto the SET island, and are 
described in more detail in appendix IE II In this manner, the two subspaces may be distinguished by measuring the 
SET current. 

Since the measurement is QND in the triplet subspace, we use the same arguments as section UlI CI to estimate the 
SET current when the DWS is a triplet. Assuming electrons tunnel between the leads and the SET island at a rate 
7, then for the configuration of lead energies described above, the current for the triplet state will be it = 67/2. 

The singlet state is approximately an eigenstate of iJnub and /Jmeas; so the same reasoning concludes that the 
singlet current should be approximately zero, is ~ 0. The measurement time for distinguishing these two currents is 
then roughly ^mcas = 2/7, just as in section llll CI 

As mentioned above, the singlet state is not quite an eigenstate of the system or measurement Hamiltonians, so 
there are corrections to the latter part of this argument. The dynamics mix the singly occupied state \S) and the 
doubly occupied state \D^). Thus there is a small amplitude for the evolution to induce transitions from \S) to 
In general these transitions are strongly inhibited since there is a large energy gap ^ U to excite the DWS to the 
doubly occupied state. Appendix IE 21 shows that in the steady state, the probability for the DWS to be in a singlet 
state is very close to unity, (S'|p'**|5) — {S\pff + pf^\S) = 1 — (^ + A)^/C7^. This means that the measurement on the 
singlet subspace is indeed almost QND, since the singlet is not changed greatly during measurement. Associated with 
the infrequent fluctuations between 15) and |-D+) is a small current. To estimate an upper bound on the current, is 
that could flow through the SET when the DWS is in the singlet state, we compute the rate at which electrons cycle 
on and off the SET island. In appendix IE 31 we show 15 < 2e7 — W — — <C it- This shows that is and ir are very 



different, and the measurement is indeed close to QND. 

C. Eff'ect of asymmetry 

Asymmetry, e 7^ 0, couples the states |I?+) and |-D-), evident in the form of iJmcas, in Eq. (|24|) . We find the steady 
state probability for the DWS to be in a singlet, (^IpIS*) = 1 — {S + A)^/(2[/^). We also compute the rate at which 
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electrons cycle on and off the SET island, described in appendices IE 31 and IE 41 which leads to an upper bound on the 

current given hy is < e{'fi + 7^)-^^ —. Both the steady-state singlet probability, and the upper bound on the 

current for the asymmetric case are similar to the results for the symmetric case derived in appendices IE 21 and IE 31 
indicating the measurement is rather insensitive to SET asymmetry. For larger values of e, when rj^lrjy'p ^ (^I^)j 
the DWS is driven into the doubly occupied subspace. which is the basis of an alternative singlet-triplet measurement 
scheme [T^ . 



D. Summary 



In principle the scheme outlined above enables measurement in the singlet-triplet basis. We have shown that it is 
possible to measure distinguishable currents through the SET depending on the state of the DWS. The energy scales 
for the singlet-triplet measurement are smaller than those for the single-particle system, by a factor of {S + A)/U, 
and this requires the lead temperatures to be smaller by a similar factor. We have also established that asymmetry 
in the fabrication of the device is less of a problem for this measurement scheme than for the single particle scheme. 



V. ESTIMATION OF EXPERIMENTAL PARAMETER VALUES 



We now estimate the required parameters for various experimentally realisable systems. Firstly we have assumed 
that 6, A U . We also require that the temperature be smaller than the SET energy shift induced by the DWS. 



Thus, for the single particle case, we require that ksT <g; 5, and Fig. 5(a) indicates that for the two-particle system, 
the temperature must be smaller than the splitting between 102 and wa, i.e. /cbT ^ 5{5 -\- 2A)/U . This is obviously a 
tight constraint on the temperature of the system. Furthermore, as discussed in section El we require the line- width 
of the SET island state must be smaller than the energy level splittings, i.e. 7 ^ S for the single-particle case and 
7 <C S{S + 2A)/U for the two-particle case. Finally, cotunneling will contribute a background current due to tunneling 
via virtual population of the SET island. 

In order to estimate the various parameters introduced for this problem, we need to estimate the overlap {L\R). 
We approximate the localised states as s-orbitals bound to each site, so that (r|_L) = ^^g-pk-ri,! y where /j, 
is the inverse Bohr radius. The integral {L\R) may be performed in prolate ellipsoidal coordinates [lj| to give 
{L\R) = {1 + fj,d + fPd? /3)e^^'^, where d = Ir^ — r/fj is the separation between the double well minima. 

We can estimate A, the tunnelling rate between localised states, for the case of an s-orbital bound to a donor atom, 
and we will use this estimate for the case of electrostatically defined gates as well. Following a similar argument to 
the derivation of Eq. we can show that A ^ 2{{R\H\L) - {L\R){L\H\L)), where H ^ pV(2me) -f VL(r) + VR(r) 

is the time-independent, single-particle Hamiltonian. Assuming Vl.r ~ q/\rL^ii\, we may again evaluate the integrals 
{R\H\L) and (L\H\L) in prolate ellipsoidal coordinates |Q| to find that the single particle tunnelling rate is given to 
reasonable approximation by A « |g/i^(ie^'^^. 

For many materials, e.g. Si or GaAs, q ~ 0.1 eV nm. For a donor atom system in Si, where |rs| ~ 10 nm 
is a reasonable estimate for the height of the SET island above the donor system, and d « 30 nm and a Bohr 
radius of = 3 nm, giving (L\R) = 2 x 10~^, and ignoring the anisotropic effective mass of Si [THl H^. we have 
6 w q{L\R) /\rs\ ~ 20 ^cV. Similarly A « 15 /icV. These figures could be increased to perhaps 200 fieV using 
an external J-gate, since the overlap integral depends exponentially on the J-gate potential. Finally we estimate 
U ^ 10 meV [i^llfl for donor impurity systems. In an electrostatically defined system such as GaAs dots, reasonable 
estimates for the various parameters are A « 150 /xeV, S « 100 ^eV and U 1 mev [TolllGj. 

For example, suppose [/ ~ 1 meV and A « 5 ~ 100 /ieV, so A/U ~ 0.1 then for the two particle case we have 
ksT <C 100 /ieV, i.e. T <C 1 K. Therefore, it is conceivable that the singlet-triplet measurement could be done at 0.3 
K in electrostatically defined dots, which is an accessible electronic temperature. In double donor systems, such as 
the Kane proposal A w (5 ^ 100 jicV is still reasonable, but since U ^ 10 meV, the relevant temperature is around 
ten times smaller, which is probably too small to be practical. This problem would be resolved if a sufhciently large 
J-gate voltage could be applied to increase A and S. Assuming that 'yi = 'jr = ksT = 30 /icV ~ 10^^^ s~^, then 
It ~ 1 nA. The proposals arc likely to work at temperatures ksT ~ 5, A as well, but with faster mixing times and 
longer measurement times. The fundamental requirement for both of the proposals in this paper is that the overlap 
between the localised wavefunctions, {L\R), be as large as possible, and preferably as large as about 0.1. 

Finally, we estimate the effect of cotunneling by comparing the conductance due to resonant tunnelling processes, 
Gres) with that due to cotunneling, Gcot- For weak coupling between the leads and the SET island these quantities 
are given by 0, 0| 

_ GlGr ttHGlGr (kBTY 

Gl + Gr ie^ 0'^ 
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where Gl and Gr are the conductances of the left and right SET-lead tunnel barrier. For the sake of estimation, 
we assume that these are equal to Gl = Gn = S^e^ /h with ^ <C 1 for weak coupling. The additional current due to 
cotunneling is small as long as Gcot ^ Gros, i-C when /csT ^ which is a less stringent constraint than above. 

Therefore, as long as the previously discussed conditions are met, cotunnelling is small. 

VI. CONCLUSION 

In this paper, we have presented and analysed a proposal for performing measurements in non-localised bases of 
both singly- and doubly-occupied double wells. 

The physical mechanism by which the measurements operate is to detect the small variation in electronic charge 
density near the midpoint of the DWS. Based on reasonable estimates for the system parameters the small difference 
in the Coulomb potential at an SET island due to the different charge distributions of different states of the DWS is 
in principle detectable. The detected signal is the current through the SET island. 

The main difficulties in these schemes is the precision with which the SET island must be placed at the mid-plane 
of the DWS and the required operating temperatures. For the single particle parity measurement, a misplaced SET 
island produces a measurement in the localised basis, which has been discussed in the past . 

The two-particle singlet-triplet measurement scheme is less sensitive to asymmetry in the placement of the SET, 
but requires very low temperatures to work effectively. 

The required tolerance to such misplacement is at the edge of current fabrication technology of 1 nm or less, and 
may be achievable in light of recent experiments [TH . Other constraints such as operating temperature and tunnelling 
rates are experimentally achievable. 

TMS thanks the Hackett committee, the CVCP and Fujitsu for financial support. SDB acknowledges support 
from the E.U. NANOMAGIQC project (Contract no. IST-2001-33186). HSG acknowledges financial support from 
Hewlett-Packard . 



APPENDIX A: MASTER EQUATION FOR SYMMETRIC SINGLE PARTICLE SYSTEM 

The Hamiltonian for the complete system of double well, SET and leads is given by 

Hrot = Ha^b -I- -I- -ffisland + ^^loads + Htun- (Al) 

Following the derivation of Wiseman et al. [l9| . we can write down a master equation for the reduced density 
matrix, R, for the system consisting of the double well plus the SET island. The dissipative terms are formally the 
same, where we identify c+ and Cg respectively with ci and b in their notation. The result is 

R{t) = — i[Hliuh + Hisland, R{t)] 

+ {7,(1 - fiiuo)) + 7r(l - /r(^o))P[(l - n+)c,]R{t) 

+ UfM + JrfrM}V[il - n+)cl]Rit) 

+ {7;(i - Mujo + 5)) + 7:.(i - /.(^o + 5mv[h+cs]R{t) 

+ HMu;o + S) + j'Jri^o + d)}V[n+cl]R{t), (A2) 

where /,(^) is the Fermi distribution for lead l{r), V[A]B = J[A]B - A[A]B = ABA^ - \ {A^ AB + BA'^A) and 7, = 
TrgilTifcgp and 7,' = 'Kgi\Tiy^^\'^ , where gi is the density of states in lead i, fco = \/2mLOQlti and fcg — ^I rajyii^ + y^jh. 
We take /z(wo) = h{^Q + (5) = 1 = fri^^o) and fr{too + = 0, as shown in Fig. |21 and then Eq. (|A2|) becomes 

R{t) = - i[Hiiuh + Hi,und,R{t)] + (7/ + irMnlclMt) + j'Mh+Cs]R{t) + j'iV[hlcl]R{t). (A3) 

We will assume that the SET island is classical, in the sense that its reduced density matrix has no off diagonal 
terms. This is justified since conservation of electron number between the leads and the island means that the electron 
number on the island is entangled with the electron number in the leads, and the lead degrees of freedom averaged 
over. Therefore we write the double well plus SET island system in the separable form 

i? = Po«)|0>(0|+pi«)|l)(l|, (A4) 

where po(i) is the state of the DWS with 0(1) electrons on the SET island. The reduced density matrix for the DWS 
alone is given hy p = Tvs [R] = + pi- 
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We now turn the master equation into a pair of coupled equations for po a-nd pi by computing the matrix elements 
{0\R{t)\0) and {l\R{t)\l) using Eq. 

Po{t) = - i[HHuh, Po{t)] - (7( +7r)-4[n_]po(i) +7r^["+]piW - 7M["+]Po(i), (A5) 
pi{t) = - i[i?Hub,Pi(<)] + ill + lr)J[n-]po{t) - j',A[h+]pi{t) + j'iJ[n+]poit). (A6) 

Since all the system operators in these equations are number operators, the equations are straightforward to solve. 
We note that po,i = o:±n± are fixed points of the equations, for some coefficients a± determined by rate balance. 
Taking 7^ = 7, and with respect to the basis {|+), |— )}, the solution to these unconditional equations is 



Poit) 
Piit) 



— 2-yt 

(A7) 



l±^p++(0) + l-^p++(0) e-3W2^+-(0) 
e-37*/2p^-+(0) e-27V„-(0) 

e-'<'/^p^+{0) (1 - e-27*)^--(o) + p--(o) 



(A8) 



where p^'' = {p\p\q). It is evident that the diagonal elements of the reduced density matrix for the DWS, p = pq + pi, 
are constant as required for a QND measurement, whilst the off-diagonal elements decay with two characteristic time 
scales, the longest of which is 2/7, consistent with the measurement time computed assuming Poissonian tunnelling 
statistics in section IlII CI 

APPENDIX B: APPROXIMATE EXPRESSIONS FOR e AND 5 

We now estimate e using Eq. (Q. For the purposes of this estimate, we will assume the SET island wavefunction 
is a delta function located at ~ {xs,ys, Zs} {xs is the axial position of the island, and Xs ^ for a symmetric 
arrangement) as shown in Fig.Ol i.e. |0s(r)p = 5{r — Tg), so 

e = Vss+- = q / d'^r- ^ - ' ' ' ' ^ 



|r-r^| |r-r,| 

where (r|L) = <Pl{^) = (f'i^ ^ tl), 0(r) is the localised wavefunction of a single site and tl is the 'centre-of-mass' 
of the left well. The last two terms cancel for all Xg, and if cc^ = then the first two terms cancel also, hence our 
previously stated result that e = for a symmetric configuration. For Xs 7^ we have 

{L\-^\L) = [ d^rML^^ = / d.d,dz , m{---L,y,z})\' 

|r-r.|' ' J |r-r,| J ^{x - x,)^ + [y - + {z - z^)^ 

We assume the asymmetry is small so that Xs is a small quantity, and we expand the square-root in a power series 
about = to find 

{L\-^\L)^{L\- L_^\L)+x,v + 0{xl), (B3) 

F fsi |r Ts I 

where r^^™ is the intended, symmetric location of the SET island and = / d'^r-^^^^^^^yi^. Following the same 
reasoning, we can show that 

We can estimate rj by assuming that the localised wavefunction is very tightly bound, so that |0(r)p ~ ^(r), and then 
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Estimating 6 is more difficult, and without detailed knowledge of the localised wavefunction 0(r) our estimate of it 
is somewhat less controlled than that of e. For a symmetric system we have 

5 = Vss++ - Vss— - (+|V^|+) - {-\V\-) 

= 2{R\V\L) -2{L\R){L\V\L) +0{{L\Rf). (B6) 

Using the approximation that (f> is tightly bound allows us to approximate ~ JFT^rv'^' ^® estimate an 

upper bound on {R\V\L) by considering that (/)^(r)(/)ij(r) is peaked with a maximum at the midpoint of the double 
well. Thus, the potential at the island due to the charge distribution (/)2(r)0fl(r) will be less than the potential due 
to the entire weight of this product located at the midpoint. That is {R\V\L) < l^jy^i . Therefore an estimate for the 
magnitude of S is 



APPENDIX C: DERIVATION OF MASTER EQUATION FOR NON-COMMUTING SYSTEM AND 

MEASUREMENT HAMILTONIANS 



In this appendix we derive a master equation for a device whose system Hamiltonian, iJgys, docs not commute 
with the measurement Hamiltonian, i^moas- The results from this appendix are used in Appendices IdI and IeI wherein 
iJsys = i?Hub. The total Hamiltonian for the device is taken to be 

^^Tot = Hgys + V + iJisland + -^loads + ^tun, (CI) 

where V = fis -ffmoas and i?isiand = t^o^s, i^ioads = J2k ^k{nik + firk) and i/tun = J2k Tikc\f,Cs + Trfc4j.Cs + H.C., as 
given in section ITTl 

The general method for this derivation follows several steps. 

1. Wc move to an interaction picture to transform away all the free dynamics. 

2. Using the Zassenhaus relation we factor the interaction Hamiltonian into a product of lead operators and a 
finite-dimensional operator acting on the DWS and SET. 

3. By tracing over the lead modes, we derive a Markovian master equation for the DWS and SET density matrix, 
R. In this master equation, the Fourier components of Bk(t) appear in the Lindblad terms. 

We transform to an interaction picture with respect to the Hamiltonian Hq = Hgys + V + i?isiand + ^icads, so that 
Htoi = Hq + Htun, and the interaction picture Hamiltonian is Hi{t) = e^"°*HTote^''"°* - Hq ^ e'"°* Htune~^"°'' ■ 
In order to compute Hi we first note that i?icads and -ff island commute with all other terms in Hq so e*^"* = 
g2(-ffays+^)tg«^fisiandtg«ffioadBt^ Uslug thc opcrator identities 

e^"c ^ c and ce^" = e^c, (C2) 

where n = c^c and \x,c\ =0, we find that 

Hi{t) = + T,fc4)e'^"''"""^*e*^^=-+^)*c,e-'(^-=+^)* + H.c. (C3) 

fe 

Since [Hsys, V]= hs® [Hsys, ^^mcas] 7^ the operator exponentials appearing above cannot be factorised. However we 
may simplify thc expression using the Zassenhaus operator relation pol |. which is an inverse relation to the Baker- 
CampbcU-Hausdorff formula, and it states that e^+^ = e^e^ Ylj] ^^'^^'^^ where each Cj[A,B] = {-iyCj[B,A] is a 
sum of nested commutators, each term of which has A and B appearing at least once (e.g. Ci[A, B] ~ ^\[A^ B] and 
C2[A,B] = — i([A, [A, B]] — [S, [B,^]]). Thc ] in thc index of the product indicates that the product is ordered in 
increasing order of j , since the factors in the product don't commute. For our purposes, the detailed form of Cj is 
not important. 
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Taking A = iHsyst and B = iVt = iHsHaiciist, and noting that fi'^ = hg so that Cj[iHsyst,insHmcast] = 
{ity~^'^nsCj[Hsys, H^eas] we have 

= g»-ffsyst -Q g-(jt)^ + ^C,[Hsys,//„ca.]g-»-ff„oastg-*-f/.y=tg^^ 
J J. 

= e'^=^=*e-'(-^=>-=+^— )*c„ (C4) 

where the first equahty foUows from direct substitution into the Zassenhaus relation, the second equaUty follows 
from repeated applications of Eqs. ljC2() . and the final equality follows by inverting the Zassenhaus relation. We may 
therefore write the interaction Hamiltonian as 

Hi{t) = 5^(T,fc4 +T,fcc|;je^('"'=-"°)*e'^=-*e-^(^=-+^— )*c, + H.c, 

k 

= ^(Ta.4+T,fe4)Bfe(i)c3 + H.c., (C5) 

k 

where 5fe(t) ^ p^{^>.-^o)t^^H,y.t^-^{H,,,+H^,^,)t ^^.g^ 

is an operator acting on the DWS alone. 

The state matrix W of the entire closed system including the double well, SET island and leads evolves according 
to the Schrodingcr equation in the interaction picture, taken to second order 

W{t + M) = W{t)-iM[Hi{t),W{t)]- M dt'[Hi{t),[Hi{t'),Wit')]]. (C7) 

Making the first Markov approximation, we assume that the system may at any time be written as W{t) = Ri{t) ® 
pi pr, so that each lead is always in a thermal state. Then averaging over lead degrees of freedom (•)/,r, and noting 

that (Q(^)fc) = (c|(^)fc) = 0, {ci(^r)kCril}k') = {4{r)k'^r{l)k') = {4{r)k'i{l)k') = ^ ^^'^ (4(r)k'^l{r)k') = " (Wfc), 

where //(r) is the Fermi distribution for lead l{r) and S{x) is the Dirac-delta function. Then the SET island plus DWS 
interaction picture density matrix, Ri{t), satisfies 



Rlit) = - I duJk {gi\Tlk\'^ fl{uJk) + gr\Trk\'^ friuJk)) X 

dt' { Bkit)Bl{t')csclRi{t') - Bl{t)clRi{t')cMt') 
~Bl{t')c\Ri{t')csBu{t) + Ri{t')csclBk{t')Blit)} 

diUk {gi\TlkWl ~ fliiOk))+gr\Trk\''il - /rK))) X 

dt' { Bl{t)Bk{t')clcsRi{t') - Bu{t)csRi{t')clBl{t') 

-Bk{t')csRi{t')clBl{t) + Ri{t')clc,Bl{t')Bk{t)] , (C8) 



where gi is the density of states for lead i. We further assume that the dynamics of the system is slow compared to 
tunnelling rates etc. so that we may make the replacement Ri{t') Ri{t) in the above integrals, making Eq. (|C8p 
local in time. Equation ljC8|l no longer depends on the lead degrees of freedom, and so is an equation for a finite 
dimensional system. With the aid of some further approximations, we may perform the integrations over ujk and 
which we now do. 

In order to do the integrals, we note that each term in Eq. IjCSp is finite dimensional so has a finite dimensional 
matrix representation. Further, we may write Bk(t) as a discrete Fourier decomposition 

N 

Bu{t)=Y,e^'^'^-^"^^'P^ (C9) 

m— 1 
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for some finite N and operators Pm ■ From Eq. (|C6p , the explicit form of Pm depends on the explicit form of ffgys and 
is important for the discussion of the dynamics of the system. The operators P„i for the single-particle DWS are given 
in Eq. I|D4|I . For the two-electron DWS, they are given in Eq. (|EH) and Eq. I)E7(I for symmetric and asymmetric cases, 
respectively. We now perform the integrations for one of the terms in Eq. HC8|) as an example, to show explicitly the 
approximations we make. For instance the second term of Eq. ljC8|l is 

dLOkgilTikWfiiLOk) f dt' Bl{t)clRi{t)csBk{t') 



mn 

~ Y,llflMP}ARl{t)CsPra, (CIO) 

m 

where 7; is defined in appendix 1X1 The first equality follows from substituting the Fourier decomposition of Bk{t), 
the second equality follows from evaluating the integral over t', and to make the final (approximate) equality we 
have made a rotating-wave approximation, where we take e^C"^"-'^™)* = Sm,,n- This is reasonable if the frequency 
difference a;„ — tOm (for n 7^ m) is sufficiently large, since when we come to solve the resulting differential equation 
terms containing a factor g^i'^n-i^mjt ^jjj ]-,g rotating rapidly, and so average to zero, to good approximation. This 
approximation is reasonable when 7;^^) is much smaller than the smallest energy level separations, 7;(,,) ^ cj„ — tjjn, 
for n ^ m. 

Applying these arguments to the other terms in Eq. (|C8|) results in the master equation for Ri{t) 

Rlii) = {Y,{llh{0J,n)+lrfrM)V[Plc\]Ri{t) 
m 

+ (7,(1 - flM) + 7r(l - frM)) V[PrnCs\Rl{t)}, (Cll) 

where again 'D[A\B = ABA^ — ^{A^AB + BA^A). This forms a generalisation of the results of Wiseman et al 
|l9l | to the situation where the measurement Hamiltonian {Hcb in their notation) does not commute with the free 
Hamiltonian of the system (i/o in their notation). Equation ljCll|) shows the importance of the Fourier decomposition 
of the system operator -Bfc(t) - the Fourier components of -Bfe(t), and their adjoint, form the Lindblad operators in 
the master equation, and it is through these components that the DWS interacts with the SET island. 
Returning to the Schrodinger picture, the master equation is given by 

R{t)^-l[Ho,R\ + Y,{{^ifi{uJm)+lrfr{uJm))V[Plc\]R{t) 

m 

+ (7KI - flM) + lr{l - fr{uJ,n)))V[P„,Cs]R{t)}, (C12) 



APPENDIX D: MASTER EQUATION FOR ASYMMETRIC SINGLE PARTICLE SYSTEM 

We now derive a master equation for the single-particle DWS for the case that the SET island is not placed 
symmetrically. Equation (|14|l gives the Hamiltonian for the Coulomb interaction between the SET island and DWS 

as 

V = ns(E) {6n+ + e{c+C- + cLc+) = fis ® -ffmcas- (Dl) 

When e 7^ 0, [-ffmoas, ^^Hub] 7^ 0, so we use the derivation of the master equation in appendix lUl 
Firstly, from Eq. (jHSl, 

Bk{t) = e'('^'"~'^°)*e*-f^H"''*e-'(^«"''+-^°™)*, (D2) 
4 

= ^ e*('^'=-'^'")*P,„, (D3) 

771 — 1 
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where a+ = = —sTa\'^)(~\ ^^"^ h± = n± ~ a±. Thus, the Fourier components, Pm, of are the operators 

appearing in Eq. (|D4|) associated with the four Fourier frequencies w„i G {ujo,ujo + S + A,ojo + S,u!o — A}. Here we 
have neglected terms of O(e^) or higher, since these are neghgible. 
The master equation for the SET and DWS is then 



R{t) = -liHiiub 



+ hi + ^r)V[dlc\]R{t) + 7;P[a_c,]i?(t) + -i[V[alc\]R{t). 



(D5) 



This expression agrees with Eq. (|ll|l in the hmit that e 0. We again assume the SET island docs not maintain 
coherence, as expressed in Eq. (|A4(I . and we then can solve Eq. (jD5() for p^{t) and pi{t). 

The most important quantity to derive from this master equation is the mixing time. By taking the Laplace 

-7/2, —87/2 and —27. All but the second of these poles appear 



transform of Eq. (|D5p . we find poles at 0, 2(5+a)^ 
as rates in the solution for the symmetric case, Eq. I|A8(I . The second pole is very small, and corresponds to the 
mixing rate induced by the asymmetry in the SET island placement. For times t <C tmix = g^^^a , the solution to 
the master equation is essentially the same as Eq. (jA8|l . On a time-scale t ~ imix, the diagonal elements also decay, 
so that the steady-state solution in the ordered basis {|+), |— )} is 



2(5+A) 









2(<5+A)''J 



and pi 



(A>A) 

1 - 

^ (5+A)^J 



(D6) 



where we have kept only the highest order term in e for each matrix element. This steady-state solution corresponds 
to the DWS being (almost) in its ground state, |— ), with the SET in the closed state, occupied by a single electron. 



APPENDIX E: MASTER EQUATION FOR A DOUBLY OCCUPIED DWS 

In this appendix we derive a master equation for the dynamics of the singlet-triplet measurement scheme. Initially 
we consider a symmetrically placed SET island. We first derive the Lindblad operators that appear in the master 
equation, and give a physical interpretation to the discrete spectrum, {wm}, that appears in their derivation. We then 
give a quantum trajectories analysis of the measured currents in the triplet subspaccs and singlet subspaces. Next we 
compute the degree of mixing in the singlet subspace induced by the fact that the singlet is not an cigenstate of the 
dynamics. Finally, we provide a derivation of the jump operators for an asymmetrically place SET. 



1. Lindblad operators 

The device Hamiltonian is given by Eq. I|25() so we can use the results of appendix O to derive a master equation 
for the dynamics of the DWS and SET island coupled to the leads. The crucial quantity to evaluate is Bk{t), whose 
definition is given in Eq. ljC6|l . Then identifying the Fourier components, Pm of Bk{t), as in Eq. (|C9p . provides the 
operators that appear in the Lindblad terms of the master equation, Eq. (jClip . 

We assume that U and a are relatively large energies, whilst e, S and A are relatively small. In fact, for the 
symmetric case, e = 0, and we will investigate this 'ideal' situation first. Clearly, if e = 0, then the dynamics is even 
more restricted, so that the measurement Hamiltonian, V = fis (S ^^mcas, only couples the spatially symmetric states 
\S) and (see Eq. (|24() '). so we may restrict our analysis to the 2-dimcnsional subspace spanned by the ordered 

basis {IS*), |Z3_|-)}. In this restricted basis we may decompose Bf^{t) into four Fourier components 

where for simplicity (and without loss of generality) we have set a = —ujq so that a common overall factor of e^i°-+'^o)t 
conveniently vanishes. Also, we have discarded terms of order 5^, (5 A and A^ appearing in the matrices, since these 
are small, but they are retained in exponents where they are the lowest order terms that lift the degeneracy of the 
tOm- We will refer to the operators appearing in Eq. (|E1|) as Pi, P2, Pi and P5 respectively. For completeness, 

Bu{i) = Pf ^''+^(t) -f e'^"' (P3 + Pe + P7 + Ps) , (E2) 

where P3 = |P_)(i:»_|, Pg = |T||)(r|||, P7 = \Tu){Tu\ and Pg |Pu>(^iil- This decomposition of Pfc(t) shows that 
there are eight Lindblad operators. Pi, Ps. 







1 _a+A 





^ 

5+A V 



U 



-iUt 



" U 





(El) 
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The frequencies Wm, associated with the measurement process, are given by 1^3,6,7,8 = 0, CJ4 = —uj2 ~ ^ ^u^^ ^^^d 
UJ5 = —LUi = U (where, again, we have set a = —ujq, otherwise we have an overaU offset of a + ujQ to our energy scale). 



These energies are shown relative to the lead chemical potentials, ^/ and /i^, in Fig. 5(a) This choice determines the 
coefficients of the Fermi factors in the master equation, Eq. IjClip . 

Wc interpret the energies fkUm as the change in energy of electrons tunnelling between the SET island and a lead. 
Thus, since 0^3^6.7,8 = 0, the corresponding processes, P3.6,7.8 are associated with elastic tunnelling between the lead 
and the SET. This can only induce dephasing of the DWS, since no energy is exchanged between the leads and the 
DWS. These elastic processes therefore do not induce mixing in the DWS, and are the origin of the QND projective 
nature of the measurement in the triplet subspace. 

Conversely, uji^2 < 0, so Pi^2 correspond to inelastic lead-SET tunnelling processes which gain an energy fiu!i^2- 
This additional energy in the lead is provided by the electron-pair in the DWS which loses energy. Similarly, processes 
P4,5 correspond to lead electrons losing energy as the DWS becomes excited. We therefore expect that there will be 
some measurement-induced energy relaxation associated with the measurement of a singlet state. 

We note in passing that these elastic and inelastic processes in the detector have counterparts in the measurement 
of a DWS by a point contact detector, as described in [T3l |. 



2. Master equation in singlet subspace 



The singlet and triplet subspaces arc not mixed at all by the dynamics, so we derive a master equation for the 
state of the DWS and SET in the singlet subspace, R ~ po <E) |0)(0| + pi (E) |1)(1|, using the results of appendix lO In 
particular, the P„i that appear in Eq. (jEl|) form the Lindblad operators in Eq. IjClljl and the cUm appear as arguments 
to the Fermi functions in Eq. I|C11|) . 

po = -i[HHuh,Po]+l'{J[P4]pi+J[P5]pi-A[Pl]po-A[P^]po), 

Pi = -i[HHuh + Hn,,,,,pi]--f'{A[P4]pi+A[P5]pi-J[Pl]po-J[P^]po), (E3) 

where 7' = 7; + 7,., and we have used T>[B]p — J[B]p — A[B]p. The steady-state probability for the system to be in 
the singlet state is given by (S'|p^''|S') = {S\pff + pV\S) = \ — {6 + A)^/J7^, is very close to unity. Therefore, if the 
DWS starts in a singlet state, its state does not change significantly during the measurement. This further justifies 
the assertion that the measurement is nearly QND on the singlet subspace. 

The poles of the master equation determine the measurement and relaxation rates. There are poles at 
0, — 7'(5^/[/^, — 7'/2 and —7'. The second of these corresponds to energy exchange processes generated by the operators 
Pi and P5 appearing in the Lindblad terms. There is therefore a measurement- induced mixing time tmix = -^iji- This 
mixing time is due to the fact that the singlet state is not quite an eigenstate of either ffmoas nor i/nub- The mixing 
time is very long compared to the measurement time, ~ I/7', since 6 <^ U. The mixing only induces relaxation of 
the DWS, and so it has very little effect on the singlet state, which is already very close to the ground state. We 
therefore conclude that this intrinsic mixing is negligible. 



3. SET average currents 

To analyse the evolution of the measurement more formally, we unravel the unconditional master equation Eq. (|C12|I . 
and derive the conditional dynamics of quantum trajectories. From the unravelling we can provide estimates for SET 
currents. Wc assume the system may be described by the density matrix given in Eq. (|12fl . i.e. the SET island docs 
not support coherent superpositions of and 1 electrons. We may therefore reduce the master equation given in 
Eq. ljC12|) to a pair of master equations for po and pi . 

The dynamics of the system decouple, depending on the state of the DWS. In particular, the triplet states and 
\D-) are eigenstates of the evolution operators, so we may consider the dynamics separately in each of the uncoupled, 
ID subspaces {|T||)}, {|T||)}, {jTn)}, In these subspaces the Hubbard Hamiltonian is proportional to the 

identity so the reduced master equations for each subspace are of the form 

Po - 7rJ[P]pi - 7iA[P]po (E4) 
Pi = 7iJ[P]po-7rA[P]pi, (E5) 

where P is the projector onto the subspace, e.g. for the subspace {|D_)}, P = \D_){D_\. 

These reduced master equations depend on only a single jump operator P, so the evolution between jumps may be 
written as a pure state, |'0ns(i))c, ris = or 1, governed by the non-Hermitian effective Hamiltonian, 7i„^ according 
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FIG. 4: (a) Projection of the state of the DWS onto the singlet state, \{S\^{tj))\, immediately after jump j for S = A/2 = 
IO7' = U/10. (b) Histogram of waiting times between jumps with = and (c) with — 1, for a simulation with lO'' jumps.. 



to the Schrodingcr equation [21| 



^MnAt))c = -iUnJi^nMc- (E6) 

at 



The non-Hcrmitian Hamiltonians for each subspacc are Tig = "^i = —i-^P. The solution is simply 

|^oW)c = e-^*|^o(0)> and IMt)), = e-^*|V'i(0)). The jump rate is determined by the cumulative density function 
(CDF) for the waiting time between jumps Vn,{tj < t) = 1 — c{'>Pn,{t)\''PnAt))c so Vo{tj < t) = 1 — e~'^'* and 
Vi{tj <t) =1 — Thus we have a cycle of electrons hopping onto the SET from the left lead at a rate 7; then 

hopping off to the right lead at a rate 7^. We therefore expect a current It = 6(7^^ +77^)^^ to flow through the 
SET when the DWS is in the triplet subspacc, in agreement with above. 

We now turn to the more complicated dynamics in the singlet subspacc, {\S), To simplify the analysis of this 

system, we will ignore the Lindblad terms depending on Pi and P5 in the master equation, Eq. ljE3p . This is reasonable 
since these terms are 0{S^ /U^), which is small. For the resulting approximate form of Eq. (|E3p the dynamics between 
jumps are governed by effective Hamiltonians for each SET island state, = 0, 1, with a single jump operator. 
Ho = H^iuh — i^-f2-P2^ and TCi = -ffnub + -f^meas — i-jPlPi- Since there is only a single jump operator associated with 
Wris , we unravel the master equations as non-Hermitian Schrodinger equations for the pure, conditional, unnormalised, 
between-jump state-vectors, \'ijjnAt))c = (OIV'^^ (0)), where ZY„^(i) = e~*^"=* During a jump at time tj, the 
state of the system evolves discontinuously according to \ipi{tj)) = P2\ipo{tj)) and IV'o(ij^)) = Pili'iitj)} ■ 

To derive an upper bound on is, we calculate the rate at which electrons hop on and off the SET, given that 
the DWS begins in a singlet state. The jump rate is determined by the CDF for the lifetime of the SET state 
with ris or 1 electrons, VnA^J < t) ~ 1 — c{^nA't)\i^nAi))c- This quantity depends on the state of the DWS 
immediately after the most recent jump. |?/'„^(t^)), which is not deterministic due to the stochastic nature of the 
trajectory. However, as discussed in section lE^ the steady state of the DWS, p'**, is very close to the singlet state since 
(S'|j0**|5) = 1 — ((5 + Ay/U^ [23. Therefore, for the purposes of computing the CDF, to very good approximation, 
we can make the replacement IV'ns(^j-)) The unnormalised, conditional state of the DWS betweens jump is 

then 

\Mt))c = Z^o|5)«e(*^-^'/2)t|5)+o(A/[/), 

\Mt)), = Ui\S) « e(*^- ^'^lut^'' +0{{5 + A)/U) , 



where t is measured from the previous jump, ft follows that the CDFs for the jump times are Vo{tj<t)Kil — e * 
and Vi{tj < t) sa 1 — e ^ . These CDFs show that the lifetime of an empty SET is short, tq = I/7', whilst 

the lifetime of an occupied SET is very long, ti = 'x^^s+a)^^' ^^^^ cycle time for electrons to hop on and off 
the island is approximately ti, and this provides an upper bound on the SET current when the DWS starts in the 
singlet state, is < e/ri. This is much less than i^, in agreement with the heuristic prediction that is ~ 0. 

To confirm these approximate analytical predictions, wc have performed numerical simulations of the conditional 
dynamics derived from the master equation Eq. (|E3p . which is not subject to any of the approximations made in this 



section. In Fig. 4(a) wc show the projection |(5|i/'(tj))| immediately after jump j for a sample of 200 jumps, assuming 
parameter values d = A/2 = IO7' = U/IO. In this figure it is evident that the DWS typically remains close to the 
singlet state after every jump, in agreement with the preceding analysis. 

Furthermore, we plot the distribution of waiting times between jumps for a simulation of 10000 sequential jumps, 
shown in Fig. 4(b) for an empty island (n^ = 0) and in Fig. 4(c) for an occupied island (ris = 1). For the parameters 



chosen, the figures show that the SET remains empty for a characteristic time of around 100/[/, whilst the typical 
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FIG. 5; (a) Transition energies, u>m, relative to lead Fermi levels, for a symmetrically placed SET island, where uji = — ojs = 
— U,L02 = —lua = —(5^ + 25/S.)/U and 0^3,6,7,8 = 0, and (b) for an asymmetrically placed SET island, where lji = —105 = 
-U,LU2 = -(5^ + 2(SA)/C/,aj3 = -c^4 = -2e and c^6,7,8 = 0. 



occupation time is around 28000/^7. These times are in agreement with the analytic estimates given above for which 
To = 100/C/ and n = 28000/[/. 

This analysis establishes that since the DWS remains close to the singlet state at all times, the measurement is 
close to an ideal QND measurement. 



4. i3fc(t) for asymmetric DWS 

We now consider the effect of asymmetry in the placement of the SET island, where e 7^ 0. As in appendix O we 
express Bk(t) in a Fourier decomposition restricted to the ordered basis {\S) ,\D^) ,\D^)}, which spans the singlet 
subspace. The Fourier components of this operator are crucial for deriving the master equation with which to analyse 
the system. We find 
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(E7) 



To derive this result, we have assumed that {5 + A)^/4[/ e ^ U. In the opposite limit, {S + A)^/4[/ 3> e, the effect 
of asymmetry is negligible. The energies LOm appearing in the exponents above are shown schematically in Fig. 5(b) 



As described in appendix IF II these energies correspond to the change in energy of an electron as it tunnels between 
a lead and the SFT, gaining or loosing energy as it interacts with the DWS. 
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